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Abstract 

First-principles calculations were performed to investigate the ferroelectric proper- 
ties of barium titanate and bismuth fcrritc , as well as phonon dispersion of BaTiOa 
, using density functional theory and density functional perturbation theory. Results 
show that the strong hybridization of Ti-0 and Bi-0 lead to the corresponding 
mechanisms for stabilizing the distorted structure. The spontaneous polarization 
of 59.4 iiCjcvr? and 27.6 iiCjcrr? were calculated for BiFeOs and BaTiOs re- 
spectively, using berry phase method within the modern theory of polarization. 
The stereochemical activity of Bi-6s long-pair, which was the driven mechanism 
for ferroelectricity in BiFeOs, was able to produce greater polarization than the 
Ti off-centring displacement in BaTiOs. New multiferroic perovskite type materials 
combined with these two ferroelectric instabilities were predicted to have a better 
ferromagnetic ordering in comparison with BiFeOs. 
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1 Introduction 



Barium titanate (BTO) is a classical example of ferroelectric material which 
was widely used in electronic devices as high permittivity capacitors, infrared 
detectors or transducers due to its particular characteristics [1-6]. The per- 
ovskite structure materials have been intensively investigated for half a cen- 
tury. Much experimental works have been conducted to analyze the properties 
of the perovskites type ferroelectric materials in order to push forward the ap- 
plied progress. Meanwhile numerous theoretical works have been employed to 
discover the driven mechanism of the ferroelectric properties in BTO, which 
is caused by the Ti off-centring displacement. 

First principles calculation has been used to describe the electronic structure 
of materials, and one of the first theoretical investigations of ferroelectric tran- 
sition in BTO has been done by Cohen and Krakouer in the early 1990s using 
the first principles calculations based on density functional theory (DFT) with 
local density approximation (LDA) method [7,8]. 

In recent years, magnetoelectric materials have attracted much interests again 
due to the large electric polarization and saturated magnetization measured 
in BiFeOs thin film. Previous works have shown that the magnetoelectric 
properties have been originating from the coupling between ferroelectric and 
ferromagnetic order parameters [9- 13]. BiFeOa has long been known to be fer- 
roelectric with a Curie temperature of about 1103 K and antiferromagnetic 
with a Neel temperature of 643 K. The Fe magnetic moments were coupled 
ferromagnetically in the (111) plane and antiferromagneticUy in the adjacent 
plane along [111] crystalline direction,which is called the G-type antiferromag- 
netism . It also has a rhombohedral distorted perovskite structure with space 
group R3c .The ten- atom rhombohedral unit cell can been seen arising from 
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a distortion of a reference Pm3m cubic perovskite structure, which has been 
considered as the probable paraelectric(PE) phase. The symmetry permits a 
canting of antiferromagnetic sublattice resulting in a weak ferromagnetism. 
However, there is a spiral structure in which the antiferromagnetic axis ro- 
tates through crystal with a long-wavelength period of 620A[14]. Therefore, 
the linear magnetoelectric effect was reduced with the cancellation of magne- 
tization due to this long spiral structure. This cancellation can be suppressed 
in thin film as a large saturation magnetization (~ 1 /ib per unit cell) have 
been observed in epitaxial BiFeOs thin film of 70 nm [15]. 

We have previously investigated the coupling of the weak ferromagnetism and 
ferroelectricity, and the ferroelectric mechanism in bulk BiFeOs , which is 
achieved by the Bi-6s stereochemical activity of long-pair. The ferroelectric 
mechanisms for these two perovskite type materials are different , which were 
mainly achieved by the different A and B cite in the perovskite structures. 
Although different works have been performed to specify the electronic struc- 
tures and ferroelectric behaviors of ABO3 perovskites materials, unfortunately 
fewer literatures can be found on the comparison of the ferroelectric mecha- 
nisms in these two materials with different ferroelectric mechanisms. In this 
paper , we mainly focused on the calculations of orbital-resolved density of 
states (DOS) ,the band structure, and the charge density of these two materi- 
als, in order to gain an clear understanding on these two driving mechanisms in 
perovskite type ferroelectrics. Moreover phonon dispersion of BaTiOa was re- 
ported here to give a further understanding on the destabilizing modes , which 
related deeply with the ferroelectric behaviors. We hope that the macroscopic 
magnetization of BiFeOa can be improved by doping with Ti. Doping of Ti in 
BiFeOs would generate the coexistence of excellent ferromagnetic and ferro- 
electric properties in the same phase. 
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The remainder of this article is structured as follows:In section II, we presented 
the computational details of our calculations. We provided the calculated re- 
sults and discussions in section III. In section IV , the conclusion based on our 
calculation were given. 

2 Computational details 

2.1 Methodology 

Calculations in this work have been done using the Quantum-ESPRESSO 
package[16], which is based on the density functional theory(DFT) and density 
functional perturbation theory(DFPT) using the plane-wave pseudopotential 
formalism. We used our self-interaction-corrected ultrasoft pseudopotential 
implementation with the Perdew Burke Ernzerhof (PBE) exchange correlation 
functional, as the common local density approximation failed to obtain a band 
gap in the transitional metal oxides. In our DFT computations, plane-wave 
basis set with kinetic energy cut-off of 40 Ry was employed, and it has shown 
that the results are well converged at this cutoff. Bi 5d ,6s, and 6p electrons, Fe 
3s, 3p, and 3d electrons, Ba 5p, and 6s electrons, Ti 3s, 3p, and 3d electrons, and 
O 2s and O 2p electrons have been treated as valence states. We used up to 
6x6x6 grids of special k points in total energy calculation. In the local 
spin density approximation plus U(LSDA-I-U) framework, the strong Coulomb 
repulsion was treated by adding a Hubbard-like term to the effective potential, 
leading to an improved description of the interaction in the transitional metal 
oxides. The Hubbard parameter U in the range U = eV - 8 eV were used 
in this work. The band gap of BiFeOs was opened significantly when U was 
equal to 7 eV , while the insulated property has not been improved as U was 
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increasing further. We used this value in the following sections.The phonon 
dispersion curves were computed within a variational formulation of density 
functional perturbation theory (DFPT) [17]. 

2.2 Structure optimization 

For BaTiOa, The lattice parameters , ranging from 3.69 to 4.29 A , was eval- 
uated with respect to the lowest total energy , which covered the average 
value(a=3.9904A, c=4.0689 A) calculated from our experimental samples, 
BaTiOs powder. Meanwhile the tetragonality maintained unchanged under 
PAmm symmetry ,as tetragonal phase is known experimentally to be stable 
at room temperature. The calculated lattice constant value was plotted with 
respect to energy and fitted to the Murnaghan equation of states , permitting 
to give the bulk modulus[18].The atomic positions were relaxed afterward by 
minimization of the Hellman-Feynman forces within a convergency threshold 
of 10~^ Ry/Bohr , using Broyden-Fletcher-Goldfarb-Shanno(BFGS) variable 
metric method, as cell shape and volume were fixed. For BiFeOa ,the lattice 
parameters and atomic positions were taken from the corresponding value of 
BiFeOs in Ref.[15]. There were smaller deviations from the original values 
after relaxation . 

2.3 Spontaneous polarization 

In the modern theory of polarization approach[19-21], the total polarization 
P for given crystalline geometry can be calculated as the sum of ionic and 
electronic contributions. The ionic contribution is calculated by summing the 
product of the position of each ion in the unit cell with the nominal charge 
of its rigid core. The electronic contribution is determined by evaluating the 
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phase of the product of overlaps between cell-periodic Bolch functions along 
a densely-sampled string of neighboring points in k-space. 

3 Results and discussion 

For BaTiOa, The lattice parameters covering the exprimental value(a=3.990A, 
c=4.069 A) was evaluated with respect to the lowest total energy . The opti- 
mization of volume was performed with lattice constant varying from 3.69 to 
4.29 A , and the volume plotted versus energy was shown in Fig.l. The cal- 
culated value was fitted to the Murnaghan equation of states, permitting to 
evaluate the bulk modulus. The normal parabolic curve gave the optimal value 
of 4.013A ,a little bit larger compared with the corresponding experimental 
value of 3.990A for BaTiOg. 

The calculated band structure of BaTiOa was shown in Fig. 2 so as to un- 
derstand the role of ferroelectric distortions on the electronic structures and 
bonding behaviors. The band gap was 1.92 eV within our approach, and this 
was in good agreement with the fact that the band gap was underestimated 
by a factor of 2 within the DFT-LDA method. The Band structure was sep- 
arated into three portions in the valence band. The narrow band positioned 
at -17eV was derived from 0-2s states, and the one positioned at -lOeV was 
originated from the Ba-5p states. There is a manifold in the valence band 
near Fermi level which was attributed to the hybridization of Ti-3d and 0-2p 
states. Moreover the conduction band near Fermi lever have a strong Ti-3d 
characteristic, indicating the configuration . There was fewer Ba states in 
these two regions ,which can be confirmed by the orbital-resolved density of 
states (DOS) of Ba 6p , Ti 3d , and O 2p shown in Fig. 3. In conduction band, 
the Ti-3d states distributing mainly in the energy region from 1.8 ev to 7.8 
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eV, overlapped deeply with the 0-2p states in the same energy region with 
slightly Ba 6s involved. The strong Ti-0 bonding under the hgand field can be 
observed from the charge density along (110) plane containing the Ti and O 
atoms shown in Fig.4. Along the displaced direction, Ti was deeply involved 
in the bonding with O to produce the covalent bonding ,which was sensitive 
to the ferroelectric displacements. The strong covalent bonding characteristic 
favored the Ti off-centering displacement along the polarized direction, which 
clearly confirmed the destabihzing ferroelectric modes, while the shghtly Ba 
involving feature tend to produce the ionic characteristic of Ba-0 bonding, 
which can be seen from the charge density along (001) plane shown in Fig.5. 
It is clear that the long-distance ionic Ba-0 bonding is slightly participating 
in the ferroelectric distortions. 

In order to gain more insight into the role of ferroelectric distortions on elec- 
tronic structure and bonding behavior in BiFeOa, a classic multiferroic with 
rhombohedral perovskite structure under R3c space group, in which the weak 
ferromagnetism and ferroelectricity coexist. We reported the band structure of 
BiFeOa in Fig. 6. The excellent magnetoelectric properties in this fascinating 
materials have attracted more interests in recent years in terms of the excel- 
lent polarization measured in epitaxial grown thin film[15]. The ferroelectricity 
is mainly attributed to the Bi-6s activity of stereochemical long-pair, which 
can only occur if the cation ionic site does not have inversion symmetry. The 
valence bands in the vicinity of the Fermi energy was originated from hy- 
bridization of Bi-6p electrons and 0-2p electrons, unhke the case in BaTiOs, 
which was caused mainly by the B cite transitional metal ion hybridizing with 
O. Moreover, the strong overlapping bands in this energy region definitely 
confirmed the long-pair feature of Bi-6s electrons, which can create greater 
ferroelectric distortions in the perovskite materials in comparison with that 
driven by the B-site off-centring in the BaTiOs. The orbital-resolved DOS of 
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Bi-6p states, Fe-3d states, and 0-2p states were illustrated in Fig.7. Fe-3d 
states were pushed away from the mixing energy region of Bi-0 in this rhom- 
bohedral structure, with fewer states taking part in the overlapping. That is 
equivalent to saying that the Fe directly connecting ferroelectricity and ferro- 
magnetism were excluded from the role of ferroelectric distortions in BiFeOs, 
which preclude the possibility of direct linear magnetoelectric effect. 

The spontaneous polarization was evaluated using the berry phase method 
based on modern theory of polarization. Born effective charges(BECs) were 
estimated by computing the Cartesian components of the polarization with 
respect to the atomic displacements, ie.: 




where A/z^^ is the displacement of ion j in Cartesian direction Zj^^ is its 
BECs tensor, and fl is the unit cell volume. The P is calculated from berry 
phase method by considering a specific structural pathway parametrized by the 
change in atomic displacement connecting a centrosymmetric reference struc- 
ture and the displaced structure. The displacements along [001] for BaTiOa 
and [111] for BiFeOa direction were chosen smaller (2% lattice constant)enough 
to ensure the validity of the linear treatment in Eq.(l). The BECs was also 
calculated based on the density functional perturbation theory[17] as follows: 

The general quadratic expression of the energy as a function of the phonon 
optic coordinates u and the electrical degrees of freedom(i.e., the field, E)is 

E{u, E) = ^Muy - ^e^E'' - eZ*n ■ E, (2) 

Where M is the nuclear reduced mass, Q, is the volume of the unit cell, e^o the 
electronic dielectric constant of the crystal, Z* is the born effective charge of 
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the ions. The electrical induction D is defined as: 

Att dE An 

Born effective charge tensor of the sth ion is the partial derivative of the 
macroscopic polarization with respect to a periodic displacement of all the 
ions of the s species at zero macroscopic electric field: 

The BECs based on these two different schemes were listed in table I. They 
yield the results in the same magnitude. It can be found that the Ti , OT(on 
the top of Ti),and Bi are all having anomalous high value in comparison 
with the nominal charges. These results agree well with the fact that the Ti 
off-centring ferroelectric distortions in BaTiOs and stereochemical activity of 
Bi-6s long-pair in BiFeOs lead to the ferroelectric instabihty. The large Ot also 
refiect that the destabihzing ferroelectric mode of Ti is along [001] direction 
in the tetragonal BaTi03,which is moving together with Ot , while Op(in the 
Ti02 termination) were slightly involved. The spontaneous polarization were 
59.4 ixC/cm^ and 27.6 ixC/cm^ for BiFeOs and BaTiOa respectively, which 
are in good agreement with the previous calculated values. The relative large 
polarization of BiFeOs indicate that the stereochemical long-pair of Bi 6s tend 
to yield a greater ferroelectric distortions than the Ti off-centring in BaTiOa. 

The calculated phonon dispersion curves along the high symmetry lines of the 
tetragonal Brillouin zone are shown in Fig. 8. The unstable modes, which de- 
termine the nature of the phase transitions and the dielectric and piezoelectric 
responses of the compounds, have imaginary frequencies. Their dispersion is 
shown below the zero-frequency line. The most unstable mode is at F, and this 
mode, dominated by the Ti displacement against the oxygens , is the one that 
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freezes in to give the ferroelectric phases. However, the instabihty is not only 
restricted to the F points. Branches of Ti-dominated unstable modes extend 
over much of the Brillouin zone. The dispersions of the transverse optic mode 
mode are flat from X to M points, and stiffened rapid towards R point in the 
Brillouin zone. This point reflect that the instability was confined to three 
quasi-two-dimensional slabs of reciprocal space, which were intersecting at F. 
This indicate a chain-like unstable localized distortions for the Ti displace- 
ments in real space, and agree well with the charge density of BaTiOs along 
(110) plane shown in Fig.4. 

The macroscopic magnetization of BiFeOa were canceled by the antiferromag- 
netic G-type structure and the spiral structure ,in which the antifcrromag- 
netic axis rotates through the crystal in a long-wavelength period of 620A. 
This drawback can be overcome partly by doping with Ti. Therefore, the 
new multiferroic combined with these two mechanisms of ferroelectricity(Bi- 
6s long-pair and Ti off-centring displacement) were expected to have excellent 
ferroelectric and ferromagnetic properties in comparison with BiFeOa. 

4 Conclusion 

Different chemical mechanisms of ferroelectric distortions in BaTiOs and BiFeOs 
have been analyzed using density functional theory and density functional per- 
turbation theory. The stereochemical activity of Bi-6s long-pair, which was the 
driven mechanism for ferroelectricity in BiFeOs, was able to produce greater 
polarization than the Ti off-centring displacement in BaTiOa. The chain-like 
Ti-0 displacement, which was determined by the strong hybridization of Ti- 
Ot, was confirmed again by the phonon dispersion curves along the high 
symmetry points in the reciprocal space. The different A and B-cites ele- 
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ments ,and structures lead to these two mechanisms of ferroelectricity. The 
Ti-doping BiFeOa was expected to be a potential multiferroic having excellent 
ferromagnetic and ferroelectric characteristics. 

This work was supported by the Aeronautical Science foundation of China 
(Grant Nos. 2003ZG51069). 
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Figure captions: 

Figure 1: Volume vs. total energy per formula unit for BaTiOs . 

Figure 2: Electronic band dispersion for BaTiOs. The Fermi level was set to 
zero. 

Figure 3: Orbital-resolved DOS for Ba,Ti, and O in BaTi03.(a),(b),and (c) 
show the DOS for Ba-5p,Ba-6s,Ti-3d, and 0-2p states respectively. Spin- up 
states are shown in the upper portions and spin-down states in the lower 
portions in all panels. The Fermi level was set to zero. 

Figure 4: Charge density (in arbitrary units) in (110) plane for BaTiOa. 

Figure 5: Charge density (in arbitrary units) in (001) plane for BaTiOa. 

Figure 6: Electronic band dispersion for BiFeOs. The Fermi level was set to 
zero. 

Figure 7: Orbital-resolved DOS for Bi,Fe, and O, and total DOS in 
BiFe03.(a),(b),(c),and (d) show Bi-6p states, Fe-3d states, 0-2p states, and 
total DOS respectively. Spin-up states are shown in the upper portions and 
spin-down states in the lower portions in all panels. The Fermi level was set 
to zero. 

Figure 8: Phonon dispersion curves for BaTiOs along the high symmetry 
lines in reciprocal space. 
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Table 1 

The polarization(P) and BECs(Z*) for BiFeOa along [111] and BaTiOs along [001]. 
The subscript AP and DFPT refer to values calculated from the approximation and 
DFPT approach respectively. 





Bi 


Fe 





Ba 


Ti 


Or 


Op 


^ AP 


4.41 


3.32 


-2.58 


2.48 


6.85 


-5.38 


-1.98 


'^*DFPT 


5.04 


3.16 


-2.20 


2.51 


7.10 


-5.56 


-2.08 


P{fi C/cm2) 




59.4 






27.6 
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This figure "Fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0711.3942vl 



This figure "Fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/0711.3942vl 
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